00 

o 

O 
(N 

>>■ 



in 



O 

o 



in 
o 

00 

o 



Fermi's golden rule and exponential decay as a 

RG fixed point 

Edwin Langmann* and Goran Lindblad^ 

Department of Theoretical Physics 
Royal Institute of Technology 
SE-106 91 Stockholm, Sweden 

May 5, 2008 



Abstract 



We discuss the decay of unstable states into a quasicontinuum using models of 
the effective Hamiltonian type. The goal is to show that exponential decay and 
the golden rule are exact in a suitable scaling limit, and that there is an associated 
renormalization group (RG) with these properties as a fixed point. The method is 
inspired by a limit theorem for infinitely divisible distributions in probability theory, 
^ ' where there is a RG with a Cauchy distribution, i.e. a Lorentz line shape, as a 

, fixed point. Our method of solving for the spectrum is well known; it does not 

involve a perturbation expansion in the interaction, and needs no assumption of a 
weak interaction. We use random matrices for the interaction, and show that the 
ensemble fluctuations vanish in the scaling limit. Thus the limit is the same for every 
model in the ensemble with probability one. 



1 Introduction 



^ ■ The standard textbook derivation of Fermi's golden rule starts from a perturbation ex- 

pansion of the unitary evolution, keeping the lowest nontrivial order, and then sums over 
a dense set of final states, to get a decay or reaction rate for an unstable state 

r = 2TTpv^ (1) 

Here v represents the (average) transition matrix element and p is the density of final 
states. (We will use the convention h = 1 throughout.) The formula is essentially 
contained in Dirac [8], see also [23]. If the rate F is a constant there results an exponential 
decay of the occupation number p{t) = exp(— Ft)p(O). The corresponding quantum 
amplitudes are the Fourier transforms of a Lorentz line shape function, see equation ([5|). 

It has been known for a long time that there are models where the golden rule and 
the exponential decay can be obtained without a perturbation expansion. It is our goal 
to show that such results hold for a large ensemble of models, and are exact in a suitable 
limit which leaves F invariant. This paper treats only the mathematical aspects of the 
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models, mainly using known methods. It does not attempt to justify their physical 
relevance, but there are numerous applications of such models to lineshape problems, 
quantum chemistry, and many other fields; we will give a few references later. 

The approach we use here was inspired by an important theorem of probability theory. 
The Lorentz line shape function ([5]) is a Cauchy probability density (PD), an infinitely 
divisible stable distribution. It is the distribution of a properly scaled limit of an infinite 
sum of independent identically distributed random variables, see Feller's book [11]. For 
any initial distribution the limit depends only on two real parameters, (F, a) in (l5|). This 
is similar to the role of the normal distributions in the central limit theorem, but the 
limit is adapted to PDs with long tails. 

A physicist would call the Cauchy PD a fixed point for a renormalization group (RG) 
of transformations on the space of PDs, see Sections [2] and[3j Our thesis is that this is also 
a good way to approach our problem. We will define an ensemble of Hamiltonians where 
the spectral density defining the decay is a closed expression (j23p in the Hamiltonian. 
There is no need for a perturbation expansion! In a simple scaling limit this spectral 
density will converge to one of the Cauchy type, a limit that can also be interpreted as 
the fixed point of a RG. 

Many attempts at a rigorous derivation of rate equations (and quantum Markov 
master equations) use the van Hove limit, where there is a scaling of the interaction 
matrix elements (v) and the time (t) as follows 

V ^ Xv, X~h (2) 

where A ^ [5,6,21,22]. In this way the dimensionless quantity Ft is invariant, but 
after the scaling the decay is on the rescaled (slow) time. Before taking such a limit it 
is necessary to let the density of states (DOS) p be infinite, otherwise there would be no 
relaxation in the limit A — > 0. Here we will instead use the scaling 

V — > Xv, p X~^p (3) 

where A — > 0, F is invariant, and t does not scale. The name continuum, or statistical, 
limit is often used, but note that here all eigenstates are normalized and the spectrum 
is always discrete (a quasicontinuum) . The dimensionless number 

Nr:=pr (4) 

scales as N-p A~^A'^p, and has a natural interpretation as the number of states under 
the resonance width. It is also the number of matrix elements effectively involved in 
the transition, a natural large parameter for the limit ([3]). We could replace the small 
parameter A by A^'p . 

We will show that ([3]) is the relevant scaling both for the limit theorem for the Cauchy 
PD, and for the quantum models introduced in Section [H and that the scaling can be 
interpreted as the iteration of a RG. In the scaling limit the models will have exponential 
relaxation, and the rate will be given by the golden rule ([T]) or a matrix generalization 
of this formula. 

An outline of the contents is as follows. 

We first discuss the probabilistic properties of the Cauchy PD in Sections [2] and [3l In 
Section S] we introduce the models of the the effective Hamiltonian type and the spectral 
density for the relaxation is found, without any kind of expansion or approximation. The 
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models are defined as an ensemble of random matrices; in the simplest case the ensemble 
is defined by just two variable (non-random) parameters p and v. 

The simplest case, a single level decaying into a quasicontinuum, is covered in Sec- 
tion m For mathematical simplicity we perform the calculations for the case where the 
unperturbed spectrum lies in a finite energy interval A with energy range /S.E. We can 
easily calculate the ensemble averages and show that they have the desired properties 
in the scaling limit. When we let A — > M the A = limit of the spectral density is a 
Lorentzian with relaxation rate V defined by ([T|). 

In Section [6] the variance of the fluctuations around the ensemble averages are esti- 
mated and shown to converge to zero in the scaling limit. Using the Chebyshev inequality 
we then conclude that each model in the ensemble will, with probability one, have an 
exponential relaxation with rate T in this limit. We also get a measure of how good this 
is as an approximation for finite parameter values; this enhances the physical relevance 
of the models. In particular, if A^p ^ 1 and T <C AE', the exponential relaxation is still a 
good approximation. We also note that the statistical ensembles can be dispensed with 
and replaced by certain uniformity assumptions on the spectra and interaction matrix 
elements. 

In Section [7] the relation between the scaling ([3]) and the RG for the matrix models is 
analyzed. In Section [8] the unitary invariance of the random matrix ensemble is used to 
show that the results will hold for more general forms of the interaction. This invariance 
also allows us to show how the properties of the Cauchy distribution under convolution 
are related to those of the spectral density. In Section [9] we sketch how the models can 
describe several decaying states, and the relation to Fano's theory of lineshape is pointed 
out. Section [10] spells out the limits to exponential decay posed by a finite spectral range 
or a discrete spectrum. Finally there are some conclusions and additional remarks in 
Section [TTl 



2 Properties of the Cauchy distribution 

The Lorentz (Breit-Wigner) line shape functions form a family of Cauchy PDs on M, 
with parameters o G M and F > 



/ dxfcixX) = 1 



Note that the mean and variance are not defined! The distribution function for a = is 

"2y' 



ry 11 
Fciy^^) = / dxfc{x,T) = - + -arctan 



(6) 



When F ^ the limit is a unit step function, Ffj = for y < 0, and = 1 for y > 0. 
Of course -F(^(— oo,F) = and Ffj{oo,T) = 1 for all F > 0. The functions ([5]) form a 
convolution semigroup 



/ dxfc{x - a^,T^)fc{y - 02 - X, F2) = fciy - - Og, F^ + Fg) 

JR 



(7) 



see [11] section 2.4. Thus, if two independent RVs and X2 have Cauchy PDs with 
parameters F^^ and Fg, then + X2 has a Cauchy PD with parameter F^^ + Fg. For 
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independent RVs {X^^, Xg, . . . , Xj^}, all with the same Cauchy PD with parameter F, 
the sum is Cauchy with parameter NT and the average is Cauchy with parameter T 

1 1 ^ 

k=l 

where ~ indicates equality in distribution. The characteristic function (CF) is 



x{t) := / dx f(j{x — a,T) ex.p{—ixt) = exp 
Jr 



-iat T\t\ 

2 ' ' 



(9) 



The convolution ([7|) corresponds to a multiplication of the CFs. 

Let the RVs in ([5]) have a common distribution F. Even when F is not Cauchy the 
scaled sum ([8]) can converge to a limit 

X^:= lim l^jv (10) 

N~>-oo iV 

with a Cauchy distribution. Sufficient conditions on F for convergence can be found 
in [11] Section 17.5, especially in the Concluding remark, involving an assumption on the 
"tails" of F. The simplest case is that the following two limits exist and are equal 

F 

lim a; [1 — F{x)] = lim xF{—x) = — (11) 
X— >oo z— >oo 27r 

where we can already identify the Cauchy parameter F. We let Xp denote the CF of F 
and define a sequence of centering constants (with a dimensional parameter 7) 



4 := 7 / sm{x /kj)F{dx), /c G Z+ (12) 
Jr 

It follows from (jlip that the following limit exists [11] 

lim [xF(t/n)e^M'' = e-l*l^/2 (13^ 

n^oo L J 

The statement for the distribution F corresponding to (|13|) is 

F^{x) := lim {*^F)[nix + (3J] = Fcix,T) (14) 



We note the similarity to the central limit theorem, but there the factor 1/A^ in (jlOp 
is replaced by l/\/iV- If F has zero mean and a finite second moment then there is a 
convergence to a normal distribution, see [11], Section 8.4. A RG version of the central 
limit theorem is described in [17]. On the other hand for an F with a finite second 
moment the limit (jl4p is a step function corresponding to a (5-function density. 



3 An example 

We now want to give an intuitive hint why the mathematical result in the previous section 
is connected to the quantum lineshape problem. As an example of a distribution which 
fulfills (jlip we pick one which is obtained from a standard lowest order perturbation 
expansion. Consider an infinite set of unperturbed quantum states {\k), k £ Z}, with 
energy eigenvalues -E^ = kuj. There is a perturbation V with nonzero matrix elements 
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(0|y|/c) = V/c 7^ 0. The perturbed state |^o) which converges to |0) as v ^ has the 
occupation numbers, to the order v'^ 



o oo 



k=l 



where v'^ must be smah enough to make Pq > 0, i.e. < Svr ^ w^. The distribution is a 
step function with steps at x = kuj 

ku!<x 



and centered, hence (1120 is zero. We find that the first hmit in (jlip reads 

2 r 

pv = — 



lim X y 775- = ^ hm X —^ = — 



.2= r 



and the second hmit is identical; consequently we know that (jl4p holds with = 0. 

Note that F * F has steps separated by lo, while (F * F)[2x] has steps separated by 
uj/2, and (*"F)(nx) in ()14p has steps separated by w/n. On the other hand we know 
from ([7]) if -F is Cauchy then this convolution and rescaling recovers F. Consequently, 
this is an example of the scaling ([3]), when we identify T with ([1]). 

As a preparation for later developments we consider a generalization to random distri- 
butions. Let V have random matrix elements (0|y|A:) = where the RVs are assumed 
to be complex- valued, normal, and independent, defined by the first two moments 

(4) = 0, (eKfc> = v'^s^, (16) 

The distribution function is now itself a random function. Consider the distribution of 
the tail function 



Tix) := 1 - Fix) = ^^Y1 



ku)>x 



Clearly the ensemble average gives back what we had before 

r 

lim X {T{x)) = — 
1^00 27r 

We can also calculate the ensemble variance. Introduce the real random vector 

Vk--=\^k\'-v' (17) 
and use the standard properties of the normal RVs ^ to calculate 

(%) = 0, iVjVk) = ^jk (18) 
Define 

AT(x) := T{x) - {Tix)) = ^ E | 

ku!>x 

Then (AT(x)) = and for large x 



(AT(x)2) 

kw>x 



1 I' dy v'^ 



yuj>x y 3lOX 
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and it follows that 

lim {/\T{xf) = lim = 

Using the Chebyshev inequality we find as x ^ oo (see [11], Chapter 5) 
Probability {IxATfx) I >5}< S^'^x^ {^Tixf) ^ > 

We conclude that with probability 1 an element F in the ensemble will satisfy pi|) . 
consequently (|13p and (|14p will hold. However, the centering constants ()12p can not be 
left out when the have the distribution defined by p6|) . 

4 Matrix models and resolvents 

In this section we introduce the matrix models of the effective Hamiltonian type. For 
most of the discussion we assume the Hilbert space to be of finite dimension. In the limit 
A ^ in ([3]) the dimension will be infinite, but again we emphasize that the eigenstates 
are all normalized. Also the limit of an unbounded discrete spectrum with finite DOS 
is interesting. We assume that these limits can be performed on the final results of the 
calculations without dealing too deeply with problems of mathematical rigor. 

We note that the formulas derived in this section are known, the main ideas going 
back at least to Feshbach [12,13,18]. The mathematical background can be traced from 
Remark 2.1 of [16]. There are also textbook treatments with physical applications, see [3], 
Complements Cj and C^^j, and [20], Chapter 21. However, we will use the results in a 
way which seems not to be standard. 

Consider a self-adjoint matrix H (Hamiltonian), with discrete spectrum {i^j^}, as- 
sumed non-degenerate for simplicity, and spectral projectors {Py} 

Define a causal resolvent (Green's function) with a regularization parameter e > 
R{z - is) := {z-ie- H)'^ = ^(z - ie - uj^y^P^ 

V 

We can pick a number of the projectors by integrating over a counterclockwise contour 
encircling the poles with real parts in a finite interval / C M: {z = -|- ie, G 1} 



y Py = — J)dzR 



2TTi 



[Z — l£) 



For small e we can approximate the contour integral by one along the real axis 

/ dx[R{x — ie) — R{x + ie)] = — [ dx Im R{x — ie) 
2m J J vr J J 

This means that we can consider the function 

- Im Rix -ie) = -Y, 7 ^ Pu (19) 

vr TT ^-^ [X — ujy)'^ + e^ 
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BjS Bj rG gularized spectral density, normalized as follows 
- / dxlmR{x -ie) = P^ = l 

and the trace of (|19p is a regularized DOS. The regularization can be written as a 
convolution (averaging) by a Cauchy PD 



Im R{x — ie) = lim / dy fc{x — y, e) Im R{y — id) 



The regularization can be justified by noting that that the duration of an observation of 
the system is bounded by a time scale 1/e, assumed much longer than the decay time. 

Now write H as a. block matrix, with projectors 11^,11^ = 1 — on the comple- 
mentary subspaces 



H 



Ha y 



B 



We know from the properties of the Schur complement [4] that 
Ra{x — ie) := n^-R(x — ie)n^ = [x — ie — Ha{x — i^)]^^ 

where, for any z £ C, we define a non-Hermitian effective Hamiltonian 
HAz) := + V{z - H^r'V^ = Hj,{z*)^ 



(20) 



(21) 



(22) 



We find for the imaginary part of (j2ip . which, with a factor l/vr, is the spectral density 
(fT9l) projected on the subspace A 



4>Aix,e) 



— Im Ra{x — is 
1 



vr 



Rjl{x — ie) e + Im H y^{x — ie) R^x — ie 



it 



(23) 



For all e > 0, x G M this is a positive definite matrix. The FT is the matrix- valued CF 
Xj^{t,e):= / dx exp{-ixt)(j)Ax,e) (24) 



(25) 



We also define the distribution function 

/X 
dy(j)Ay,e) 
-oo 

Note that the normalization implies that for every e, 



^>^(-cx),e) = 0, «>A(oo,e) 



1. 



When e — > ([25]) converges to a step function $^(x,0); the steps are a,t x = to^,, the 
spectrum of H. 

We can see in ()23p the beginning of an exponential relaxation. In fact, if we could 
assume that H^x — ie) is independent of x, then (j24p is a matrix- valued decaying 
amplitude 



XAit,^) 



exp 



(^-itReH^ - \t\e - |t|Im#A^ 



(26) 
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Our goal is to justify this simple form in the limit ([3]), where we can finally set e = 0. In 
most of the calculations below we will let the subspace A have rank 1, and then, if we 
leave out the argument x — ie, the projected spectral density simplifies to 

1 e -'rlmH a(x — ie) , , 

4>Ax,e) = = ^ '-^ 27 

^ ^ (x-ReF^(a;-ie))2 + (e + Im/7^(x-ie))2 ^ ' 

Again, if H^{x — ie) is constant we have the resonance form ([5]). 

It is important for the success of this approach that the e-averaged form of the resol- 
vent (I23p is expressed in terms of lmH^{x — ie) which is also an e-averaged expression, 
as we will see below. 



5 The simplest model 

The model (120p is particularly simple when the subspace indexed A has rank 1, and 
the spectral density has the form (j27p . This case has been solved repeatedly since a 
long time, and used for numerous applications, see for example [1,3,7, 10, 15]. Here the 
interaction V will represent an ensemble of random matrices. This is a device often used 
in applications [2,19], here it allows us to make statements valid for almost all elements 
in the ensemble. 

Now is represented by single state \tps), of energy E^, while is spanned by 
a finite set of states {\j);j = —N^, —N^ + 1, . . . ,N^}. First we choose an equidistant 
spectrum Ej = Eq + j tu^ where Eq k, E^. The calculations will show that we can cope 
with more general spectra with an average DOS = We also assume that the 

energies are restricted to a finite interval A of width AE = E_^ — E_ 

E^eA:= [E_,E^] 

Only a lower bound is essential for finite quantum systems, but here we prefer to simplify 
the mathematics by having a finite number of states in H^, = 1 + AE; later we 
can let this number approach oo. 

The interaction V couples [V's) with all the allowing this state to decay. We 
choose the matrix elements to be independent complex normal random variables 

{i^s\V\j) = {j\V\^P,)* = (28) 

defined by (jl6p . while all other matrix elements are zero. Thus, in (j20p . Hj^ = E^, V is 
a row vector with components S^j, and is a diagonal matrix with elements Ej. The 
spectral density (p7P is now 

Imi?(x - ^e) [V'.) = - , ' ^"^^'"t^^i' 2 (29) 
and in the limit e — > the FT (j24p is the amplitude for staying in the state 
(V'Jexp(-im)lV',) = ^ |(^/;Ja;^)|2exp(-itu;^) 

V 

For given S^, and V the equation system to solve for the eigenvalues oj^ and eigen- 
vectors \ljJ^) is [1] 

(i?,-^)(VJu;) + 5;e,(i|u;) = 

i 

(i?,-a;)(fc|^)+e^(V'J^) = 
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We know that each eigenvalue is located between two unperturbed {V = 0) eigenval- 
ues, but we do not need the exact values. Instead we can estimate the terms in the RHS 
of (|27p . Decompose (j22p into real and imaginary parts, for insertion in ([2j 



^ (x-E.)2 + e^ 
We first deal with the imaginary part 



lm.Hj^{x — ie) = £ 



1^/ 



(30) 



(31) 



A simpler case is solved e.g. in [1], where the interaction is non-random, = v, and 
the spectrum of is unbounded, i.e. = ±00. When to^ <C e the infinite sum is 
approximated by an integral which is independent of the argument in 
,2 



ev 



dy 



1 



y2 + £2 



r 
2 



(32) 



where we have put p = in ([T]). For our model, first apply the ensemble average over 

ippr( 

c.\2 , (33) 



the RVs ^j-, then, again, approximate by an integral 



Imif 



ev 



J(x, e) 



71" -a 



(x-EA^ + e^ 



dy 



y2 + 



arctan 



arctan 



It is easiest to see the behavior of the function J by a computer calculation. For e <C AE 
it is close to the CF x(A,x) for the interval A, with deviations e-near the end points. 
When e ^ then J(x, e) — > x(^) ^) and 



r r 

Im i/^(x - ie) ) J (x, e) ^ ^ x(A, x) - 



The same method applies to the real part 

,.2V- (X-^ ^.) 

^ (x - £;„■)' 

3 



ReH^ 



X - E-Y + £2 



+ rif(x,e) 



(34) 



(35) 



K{x,e) 



1 

2^ 



y 



y2 + £2 



47r 



In 



(E+ - X)2 + £2 



X 2 + £2 



The second term in (|35p represents an x-dependent level shift; it has a well-defined limit 
when £ ^ 0. It has the eff'ect of changing the DOS slightly, on the order of F/AS. Near 
the center x = {E^ + E_ )/2 of A and for large /S.E 



K{x,Q) = i-ln 



AE - 2(x - x) 
A£; + 2(x-x) 



2 X — X 



vr AE 



There is no unique way of taking the limit AE 
bounded, it holds for all fixed £ > that 



lim K(x,e) 

AE^oo 







(36) 

00, but as long as x and x stay 

(37) 
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Summing up the results so far: if we use the ensemble averages derived above for the 
terms in (I27p . and replacing the sums with integrals, we find the following approximation 
which we expect to be good for uj^ <C e (cf. [3], Complement Cjjj) 



1 



2e + TJ{x,e) 



27r [x-E,- TK{x, e)]^ + [2e + TJ{x, e)]74 



(38) 



6 Estimating the ensemble fluctuations 

We want to justify the calculations in the previous section by showing that the ensemble 
fluctuations and other error terms vanish in the limit ([3]), thus that Im Hj^ and HeH^ 
are given by the last terms of (j33p and (|35p . and that (j38p becomes exact. We will 
concentrate on Imi?^, since the real part is very similar. 
The deviation from the ensemble mean is given by 



Almif^ := Imi^^ - (ImH 



3 - ^' 

where rj is the RV ()17p . with mean and variance given by (jlSp. The variance of (j3ip is 



2 4 

e V 



For an order of magnitude estimate we again replace the sum by an integral, and extend 
the integral to M. Calculations similar to those performed above give 



£7^ (imH, 



vr V 



2euj 



B 



2 

2e 47r 



(39) 



For the relative size of the fluctuations we can take the square root of this expression 
over ()32p . estimating the importance of the correction due to these fluctuations by the 
dimensionless parameter 



K :- 



B 



1 



27re yj2-KN^ 



(40) 



where := p^£. From ^ follows that it scales as k ^ \k and goes to zero. Similar 
calculations show that k also measures the fluctuations in Re-ff^. We again use the 
Chebyshev inequality for an upper bound on the probability of having a deviation from 
the ensemble average [11] 



Probability{|AIm^^| > 5} < 



1 



2ttN^ 



T 
25 



n 2 



Clearly, this quantity scales as A^, and the scaling limit will give a value for Imi?^ and 
Re Hj^ equal to the ensemble average. 

We also have to deal with the error involved in replacing the sum over the spectrum 
by an integral, and in the assumptions on the level spacing. If the statistics for the level 
spacings form a Poisson process with parameter p^, then the expected number of levels 
in an interval 5E and the variance are equal 



SE/ 
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Using the random version of the sequence {Ej} in a sum hke (j30p . we get the relative 
size of the resulting fluctuations by setting 5E = e 

^ = J_ = ^ (41) 

The square root is of the same order as k (j40p . If instead we use a uniformly spaced 
spectrum, the integral approximation has a smaller error, of order (j4ip squared. We 
know that random matrix spectra are typically more uniform than the Poisson case; this 
is the feature of "spectral rigidity" due to level repulsion [19]. The correction term will 
then be in between the uniform and Poisson values. It appears reasonable to consider 
the Poisson value as a worst case for ensemble of spectra with a given density for the 
ensemble average. 

The argument so far justifies (I38p as an exact result in the limit ([3]). What happens 
when the limits e — > and A ^ M are taken after ^7 The order of taking the limits is 
not essential here. The limit e — > is straightforward in (I38p . and we get a normalized 
density </>^(x, 0). This is a modified Cauchy density in a finite spectral interval and with 
a level shift function from (|35p and ()36p . Finally, in the limit A ^ M, using (137p . we 
recover the exact Cauchy form 

lim6^{x,0) =fc{x-E„T) (42) 

In view of the fact that many applications involve models with large but finite DOS 
p^, it is also interesting to see what we can say in this case, when we let e — > while 
A^P ^1. It is then convenient to use the integrated distribution function (|25p . which 
is always finite and monotonically increasing, while (p^^ is a sum of (5-functions. Thus 
$^(x,0) is a step function, while the e-averaged version $^(x,e) is close, as measured 
by the small number (jlOD, to a Cauchy distribution ^ for 1 <^ <^ and F < AE. 
Given the properties of these two functions, the local averaging on an energy scale e 
cannot have a drastic effect, and it must hold that 

^Ax,0)^Fcix-E,,T) (43) 

is a good approximation. This statement is supported by numerical calculations, see 
Figure [TJ We can conclude that the distribution calculated from our model is close to 
Cauchy if there is a clear separation of energy scales 

« r < AS (44) 

which also implies that Nj. S> 1. 

To reach these conclusions, the properties of the spectrum of and of the matrix 
elements of V were crucial. So far we have assumed simple regular or random distribu- 
tions to do our estimates. Similar calculations are possible for more general sequences 
{Ej} and without an assumed statistical distribution. Instead we can postulate a 
uniformity for the spectrum and the matrix elements, a property which could in principle 
be verified in a concrete model. For the level spacing assume that the following estimate 
holds uniformly in x, for all e ^> uj^ and with k defined by (j40]) . 

^2 2=^ + 0i^) (45) 

TT ^ (X- E,:)^ ^ ^ ^ ^ 

7 J 
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When e ^ oo (or e AE ) the RHS must be essentially unity, and this fixes oo^. In 
the same way the uniformity of the matrix elements means that 

1^1! =i + 0{^) (46) 

where v"^ is defined through the limit e — > cxd. The factor before the sum is 2e/T. 




_10 -8 -6 -4 -2 2 4 6 8 10 



Figure 1: Illustration of the closeness of the computed distribution function (the step 
graph) to the Cauchy form given by ([6]) (the smooth graph) for a single random choice 
of the interaction V in the model of Section [5j The parameters are = 300, /S.E = 20, 
r = 1.41. Then Ny = 21, not a very large number in this context. For larger values we 
see a convergence to the Cauchy form. 



7 Scaling and renormalization 

It is easy to see how the scaling ([3]) works in the models of the type ()20p . If we disregard 
mathematical rigor the simplest form is obtained when the energy interval A = R. Then 
H in (j20p transforms in a simple way: for < A < 1 





V 






xv 













(47) 



Clearly these transformations form a semigroup, and we think of them as a RG. We can 
let and V be random or not; it is enough that they fulfill the uniformity properties 
(j45p and (j46p . which are preserved under the transformations. Of course, there is no 
proper limit for ()47p as A — > 0. For any initial choice of H of the type indicated above 
there a limit for <I>^ of the Cauchy form, and this is then the "fixed point" of the RG 
(HZl). 
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When A = M we can also discuss the relation of the RG transformation (j47p with 
that defining the Cauchy PD limit in Section [3l For large x the tail function of (I25p is 



r(x,e) 



1 _ $ 



vr 



dyy 



-2 



le] 



For simplicity, again assume that the spectrum of is {/cw^, /c G Z}. Then a straight- 
forward calculation gives 

|2 



lim Tix, e) 



fc2 



Comparing with the calculations in Section [3] we again find the tail condition, with 
probability 1 



lim xT{x,0) = — 

Hence, setting ^>^(a;,0) = F{x) we find that (jlSp and (jl4p hold. The two limiting 
procedures, that applied to the quantum models and that applied to PDs, give the same 
limits for an infinite uniform spectrum. The same result holds for a spectrum that is a 
stationary Poisson process. 

We would like to define a RG transformation also when A is a finite interval. What 
does it look like in the matrix form (j20p ? We find a scaling ([3]) with A = l/\/2, and which 
involves a doubling of the dimension of the matrices. To do this explicitly, consider a 
matrix composed of three parts, where the assumptions on Hj^ and the other components 
are as in Section [5] 



H 



Ha 



c 



H 



B 



c 








H. 



c 



(48) 



We note the matrix identity 



[VB,yc 



c\ 



H 




B 





H, 



c 



[Vi 



B 



and the corresponding decomposition of the effective Hamiltonian 



HJz) 



c) 



Thus the contributions of B and C to Imff^ and to the decay of A just add up. Clearly 
Pb-\-c ~ Pb~^ Pc- order deal as simply as possible with the energy level statistics, we 
can assume Poisson statistics for and H(j with parameters and P(j. For B C 
there will then be Poisson statistics with parameter p^^q- 

We can then let the A = 1/ \/2 transformation be represented by 

1 



Hb ~^ -f^B+B'i 



V 



V2 



[V,V'] 



(49) 



where and H^, have the same DOS p^ and independent Poisson statistics, while V 
and V both have the distribution defined in (|16p . Clearly the number of energy levels 
scales as A'^ ^ 2N. We can choose for the RG an iteration of the dimension doubling. 



hence a sequence A^ 
set e = 0. 



2 '^/^Aq —>■ 0, obtaining the limit discussed in Section [6l when we 
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8 More general models 



Here the solution of the model (f20l) in Section [5] is used to understand the properties 

of apparently more general models. We will also be able to see the relation between 

the convolution property ^ of the line shape function and the models defined here. 
Consider Hamiltonians of the form 

H^=Hq + V (50) 

where Hq is a diagonal matrix with properties like those of in (I20p . Again let 
and \k) be the eigenvalues and eigenvectors of H^, while the DOS is Pq, and the spectral 
interval A. The eigenvalues and eigenvectors of H^^^ are lv^ and The matrix elements 

{j\v\k) = {k\v\3r = i^^ 

are still assumed independent, identically distributed complex normal RVs 

{ijk) = 0> {Cjkimn) = V^^jmkn (51) 

We note the known fact that the ensemble of such random matrices V is invariant under 
all unitary transformations U, i.e. U^VU ~ V [19]. Numerical simulations indicate that 
if r := <^ AE and pgF ^ 1 then the DOS of is also pg, except near the 

endpoints of A. 

Every eigenstate of Hq will decay as a result of the interaction V. Numerical evidence 
indicates that the decay rate is near T when PqT S> 1, and that the spectral densities are 
close to those found in Section [5l If this is true, then |(A;|a;j,)p inserted in (f29l) gives a 
spectral density near the Lorentz form, for every choice of k. (This fails for E"^ near the 
end points of A, but we will argue as if A = M.) We now try to justify this picture, but 
without mathematical rigor. 

Pick any eigenvalue of Hq that is not too close to the ends of A. Call it Eq and 
calculate how the corresponding eigenstate |0) decays by transforming this problem into 
that already solved in Section [5j Make a decomposition 

where there is only one non-zero element of H^^ = E^q+^qq, Vi is defined by (V^i)ofc '■= ?ofc 
and what is left is H^^. For each ^ in the ensemble (jSip we can diagonalize by a 
unitary U which leaves the basis vector |0) invariant. The transformed Hamiltonian is 

U^H^U = H^A + U^H^bU + ViU + U^vl 

which is now of the form ()20p . where Hj^ — > is 1-D, and U^H-^b^ is diagonal with 
eigenvalues randomly distributed, but with a level repulsion which makes the DOS nearly 
uniform and equal to /9q on the average. Due to the unitary invariance of the ensemble 
(jSip . it holds that ~ V]U . Thus the solution in Section [5] will still hold here for the 
decay of the chosen state, with replaced by Pq. The fluctuation properties will also 
be the same. 

We make the following note on the invariance properties of the ensemble ()5ip . The 
following phase transformation of the interaction part 
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leaves the ensemble invariant. It is easy to check that the eigenvalues are invariant and 
that the eigenvectors are just multiplied by a phase 

Now let = {k\C,) be the components of a vector. We can allow it to be a random 
vector, but assume it independent of the ensemble defined by (|5ip . Then in 

(CK) = j;Cfe(A;K) 

k 

an average over the ensemble ()5ip . or just over all angles 9f^, gives zero, while in the 
average of the absolute square the cross terms vanish 

k 

Using the spectral density form ()29p we find for PqT ^ 1 

{mi^,)\')^Po'^ct>{u;,-E„e)\C,f 

k 

and replacing (^^ by a smooth function Ci^k) approximating by an integral 

i.e. a convolution. 

We now add another interaction term with the same statistics 

= H^ + V = Hq + V + V' 

i.e. V ~ y , but we assume that V and V' are independent RVs. By the rules of adding 
independent normal RVs with zero mean V + V' ^ \plV . The diagonalization of can 
be done in one step using the interaction V+V , or in two, first making diagonal, then 
H2, while using the unitary invariance of the ensembles. Use \^y) for the eigenvectors of 
Then expand the scalar product 

{k\nj = J2{k\^u){^M 

V 

take the absolute squared and a suitable ensemble average, as before 

V 

If we average over the angles only we need no averages in the RHS. Finally with the 
spectral density coming from diagonalizing starting from i?Q, and from diagonal- 
izing starting from we find the density for diagonalizing H^. starting from Hq as 
a convolution 

With the explicit resonance form (I38p . we find that (^^ '^^d <^2 have parameter F while 
that (\) in the LHS has parameter 2F, coming from y + y ~ \plV . Thus we recover ([7]) 
in the case F^ = F2 = F. 



15 



9 Higher dimensions and Fano lineshapes 



We return to the model (120p and allow Hj^ to have a finite dimension N^. This kind of 
model can describe several states (interacting or not), all decaying into quasicontinua. 
In such applications we must make physical assumptions on the matrix elements of V , 
and it is not always relevant to assume them all independent as in (|5ip . It is easy to 
use two different quasicontinua with different DOS, as in (f48l) . perhaps coupled to two 
subset of levels in Hj^. 

Pick one possible structure by choosing ^ to be complex, normal RVs defined by 

where 7 is a positive semidefinite matrix of dimension x N^. The imaginary part of 
()22p is then a positive semidefinite random matrix 



- ^e)]^^ = ^E (,_gj;+,2 (52) 



with ensemble average 

(Gjk) = Ijk ^ _ E,f + £2 ""PB^jk 

Comparing with the case Nj^ = 1, it is natural to define the golden rule matrix 

^jk ■= '^T^PBljk (53) 

The fiuctuations around the average is estimated as in Section [6l and again found to 
contain a factor uj^/e, just like ([39]) . Then ([3]) scales 7 as 7 ^ A^7, and we will find the 
convergence to the ensemble mean in the limit A ^ just as for iV^ = 1. The same holds 
for Rei?^. Consequently the exponential solution (I26p (including e = 0) is justified. 

When Hj^ and F commute and can be diagonalized simultaneously, each eigenstate 
decays independently and exponentially, each with a different Lorentz line shape and 
decay rate. When they do not commute, the present formalism includes models with 
Fano line shapes [7,9, 10]. Assume that we can prepare the decaying system of A^^ 
states in a definite quantum state |^). The lineshape function which will govern the non- 
exponential decay of this state is obtained from the spectral density (j23p by projecting 
it on the corresponding 1-dimensional subspace (and setting e = 0) 

f[0,x) ■.= -{d\ii^Rj,{x)\e) 

TT 

Clearly the factor i?^ in (j23p will in general have A^ poles in the upper half plane (and 
R\ the conjugate poles in the lower half plane); these poles will determine the spectral 
density, which is far from the simple Lorentz form in general. Depending on the vector 
1^) there will be a weighted interference of the contributions from the resonance poles. 
The standard Fano lineshapes can be reproduced for A^ = 2 by suitable choices of the 
parameters. When A^^ > 2 the number of different possibilities increase rapidly. 
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10 Separation of time scales 



It is well known, for general models of decay, that if the energy spectrum of the reservoir 
is bounded then the decay will be non-exponential for very short and for extremely long 
time scales; this subject was discussed in e.g. [14]. For the class of models considered 
here the Lorentz line shape and exponential decay is an exact result in the limit ^ if 
we also set the energy interval A = M. We will very briefly analyze the deviations from 
this limit caused by a finite A or a finite p^. 

First consider the case when A = M but p^ < co. When {Ef, = kuj^; k G Z}, the 
"decay" is periodic with period := ^irp^. Clearly we need to have the decay time 
defined by F to be much shorter than that. This is so if u)^ <^ F, or, equivalently A^p ^ 1. 
However, in order to find the time we need to see the deviation from an exponential 
decay to zero, another calculation is needed. Recall that $^(x,0) defined in ([25|) is a 
step function reflecting the discrete spectrum, but it still has a useful approximation 
(j43p . The decay (characteristic) function (j24p will then be almost periodic. Consider the 
ergodic limit 



lim 1 / dt\xAit)\' 



For exponentially decaying amplitudes this is zero, but for a discrete nondegenerate 
spectrum the limit picks out conjugate Fourier coefficients of Xa(^) Xa(0* reducing 
the double sum to a single (diagonal) sum which is constant in time 



1 Ll) 1 

lim - / dt\xA(t)\'^ oc = — 



When the sum is approximated by an integral we find that 
m — / 

T~ 

Thus, over long times the function |Xyi(^)P is almost periodic, fluctuating around a non- 
zero mean value determined by A^p. We can expect to see the exponential decay in this 
fluctuating background only for t satisfying 

t F < In TVp (54) 

Next we consider the effect of a finite AE while we can let p^ = oo. Then |Xyi(^)P 
has a smooth quadratic maximum at t = 0. Numerical simulations for large values of p^ 
(taking into account the level shift) show that this function is close to an exponential for 

tA^>l (55) 

but they also indicate that there is a deviation from the exponential 

|x(t)p-exp(-Ft) = 0(F/Aii;) 

which remains significant over a time interval of the order of 1/F. When F <C AE this 
deviation is very small, and the time 1/AE is very short compared to the relaxation 
time scale. Then the effect of a finite AE on the relaxation will be insignificant, and the 
exponential is a good approximation. 
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11 Conclusions 



In this paper we argue that the exponential relaxation and golden rule for model Hamil- 
tonians of the type ([20]) can be understood as fixed point properties of this class under 
a simple RG of transformations, e.g. in the form ()47p . 

The mathematical basis for the results is the simple non-perturbative form (j23p for 
the spectral density and the assumed uniformity of the spectrum of and the matrix 
elements of V, for instance in terms of the relations ()45p and ()46p . In the random matrix 
approach the results hold with probability one for each element in the ensemble. The 
method works for a single decaying state, with a simple Lorentzian lineshape, as well as 
more complex cases, including Fano lineshapes, and it does not involve an expansion in 
the coupling strength. Instead we estimated the deviations from the scaling limit when 
the scaling parameter A is small but nonzero, and found that the deviation is small if 
the number of states under the resonance is large: Nj^ = p^T ^ 1. We could also handle 
the case of a finite energy range AE for H^. Using the dimensionless quantities Nj. and 
N = p^AE (the total number of states) the condition for small deviations (I44p reads 

In Section [10] we concluded that the exponential decay is a very good approximation in 
a time interval restricted by (154p and (155p . In terms of and Nj. we found that 

Nj^/N < tr < IniVr 

On the other hand, we note that these models cannot be used to derive the most 
general irreversible evolution for open quantum systems. The CF (I26p . acting on any 
pure initial state in the subspace A, defines a pure state, with decreasing norm, for all 
t > 0. This is clearly a special case, but still useful in many applications. 

In view of the frequent use of the van Hove limit ([2]) it is important to say that ([3]) is 
not a weak coupling limit; instead it must be interpreted as a limit where the strength of 
the coupling is held constant while the reservoir is enlarged. In order to see this we first 
have to find a measure of the strength of the interaction term in the model of Section [5] 
When < oo then we can choose the square root of the positive scalar (expectation 
for a random V) 



to get a measure of dimension energy. In the scaling ([3]) the parameters F and AE are 
invariant, so is the strength of the interaction. When AE — > oo this value is inevitably 
misleading. The matrix elements connecting the decaying level to off-resonant final 
states, i.e. with |£^^ — Ej^\ 3> F, will have little influence on the decay except for very 
short or very long times. Provided F <^ AE we can then replace AE by a multiple 
of F and use the scaling invariant F itself as a useful measure of the strength of the 
interaction. 
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